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Abstract 

Background: A plant's endogenous clock (circadian clock) entrains physiological processes to light/dark and 
temperature cycles. Forward and reverse genetic approaches in Arabidopsis have revealed the mechanisms of the 
circadian clock and its components in the genome. Similar approaches have been used to characterize conserved 
clock elements in several plant species. A wild tobacco, Nicotiana attenuata has been studied extensively to 
understand responses to biotic or abiotic stress in the glasshouse and also in their native habitat. During two 
decades of field experiment, we observed several diurnal rhythmic traits of N. attenuata in nature. To expand our 
knowledge of circadian clock function into the entrainment of traits important for ecological processes, we here 
report three core clock components in N. attenuata. 

Results: Protein similarity and transcript accumulation allowed us to isolate orthologous genes of the core circadian 
clock components, LATE ELONGATED HYPOCOTYL (LHY), TIMING OF CAB EXPRESSION 1/PSEUDO-RESPONSE 
REGULATOR 1 (TOCl/PRRl), and ZEITLUPE (ZTL). Transcript accumulation of NaLHY peaked at dawn and NaTOCl 
peaked at dusk in plants grown under long day conditions. Ectopic expression of NaLIHY and NaZTL in Arabidopsis 
resulted in elongated hypocotyl and late-flowering phenotypes. Protein interactions between NaTOCl and NaZTL 
were confirmed by yeast two-hybrid assays. Finally, when NaTOCl was silenced in N. attenuata, late-flowering 
phenotypes under long day conditions were clearly observed. 

Conclusions: We identified three core circadian clock genes in N. attenuata and demonstrated the functional and 
biochemical conservation of NaLIHY, NaTOCl, and NaZTL 
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Background 

The circadian clock, entrained by our planet s 24 h rota- 
tion on its tilted axis, plays crucial roles in the 
synchronization of the performance of organisms with 
daily cycles of light and temperature, enabling organisms 
to regulate activities at the correct time of a day [1]. For 
instance, the endogenous clock in plants influences 
various biological processes including leaf movements, 
hypocotyl growth, floral transition, and abiotic and biotic 
stress resistance [2-4]. 
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The circadian rhythmicity and molecular mechanisms 
underlying the circadian clock have been investigated in 
many organisms including Drosophila melanogaster, 
Neurospora crassa, Synechoccocus elongatus, and mice 
[5-7]. In general, several interconnected transcription/ 
translation feedback loops participate to establish central 
clock oscillations [8-10]. In plants, circadian rhythmicity 
is extensively investigated in a dicotyledonous model 
plant, Arabidopsis thaliana, and a 'three-loop model' has 
been proposed [11]. 

TIMING OF CAB EXPRESSION 1/PSEUDO-RESPONSE 
REGULATOR 1 (TOCl/PRRl) and two partially redun- 
dant MYB transcription factors, CIRCADIAN CLOCK- 
ASSOCIATED 1 (CCAl) and LATE ELONGATED 
HYPOCOTYL (LHY) comprise a central oscillation loop 
[12]. CCAl and LHY repress transcript expression of 
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TOCl by binding directly to its promoter, and as shown 
recently TOCl negatively regulates the transcription of 
CCAl and LHY [13,14], establishing negative feedback 
loop [11,15]. The second loop is formed by TOCl and 
the evening complex (EC) consisting of EARLY FLOW- 
ERING 3 (ELF3), ELF4, and LUX ARRHYTHMO (LUX) 
[14]. In addition, the third loop is established by negative 
feedback between PRR7/PRR9 and LHY/CCAl [16,17]. 

Protein turnover regulation provides another layer of 
sophistication in the regulation of the circadian clock. 
Notably, the TOCl protein is regulated by proteasome- 
mediated protein degradation. The F-box protein 
ZEITLUPE (ZTL) in the E3 ubiquitin ligase SCF (SKPl- 
CULl-F-box protein) complex ubiquitinates the TOCl 
protein through direct physical interaction in a dark- 
dependent manner [18]. Homologs of ZTL are also 
involved in circadian regulation of photoperiodic flow- 
ering [19,20]. The FLAVIN-BINDING, KELCH REPEAT, 
F-BOX 1/ADAGI03 (FKF1/AD03) interacts with GI in 
a blue light-dependent manner and regulates CON- 
STANS expression by degrading a Dof transcription 
factor, CYCLING DOF FACTOR 1 [19]. In addition, it 
has been reported that ZTL, FKFl, and LOV KELCH 
PROTEIN 2 together regulate TOCl and PRR5 degrad- 
ation, contributing to the circadian oscillation [21]. 

Recent works has shown that many of these circadian 
clock components can also be found in diverse plant 
species including rice, soybean, maize, and poplar [22- 
25]. CCAl /LHY genes are widely conserved in eudicoty- 
ledonous (eudicots) and monocotyledonous (monocots) 
plants [26-29], and CCAl/LHY and TOCl feedback 
loops are thought to play a central role in the clocks 
function in these plant species. Functional homologs of 
ZTL have also been found in several plant species 
[30,31], further supporting that circadian clock compo- 
nents are fairly well-conserved in plants. 

Despite the important role of the endogenous clock 
in entraining physiological processes to environmental 
signals, how the circadian clock regulates ecological 
performance of a plant in its natural habitat is largely 
unknown. Only a few studies have shown that the en- 
dogenous clock allows plants to maximize photosyn- 
thetic capacity and reproductive success [32-34]. In 
order to expand our understanding of the clock function 
in biotic and abiotic interactions, we identified three 
core clock components (LHY, TOCl, and ZTL) in a wild 
tobacco, Nicotiana attenuata, which has been developed 
as a model system for understanding ecological perform- 
ance in native habitats, in particular the Great Basin des- 
ert in Utah. N. attenuata is a seasonal solanaceous plant, 
completing its life cycle during spring and summer to lie 
dormant in the seed bank for the many years between 
fires in its native habitat. These results provide add- 
itional evidence of the conservation of the circadian 



clock genes and set the stage for future studies to un- 
ravel the ecological relevance of the clock. 

Results 

Isolation of putative core circadian clock genes in 
A/, attenuata 

The full-length or partial sequences of three putative 
core circadian clock genes {LHY, TOCl, and ZTL) in N. 
attenuata were isolated by BLAST search against in- 
house cDNA library using the sequences of Arabidopsis 
clock genes. We checked the diurnal expression of these 
transcripts in our time series microarray database [35], 
which examined patterns of transcript accumulation in 
N. attenuata leaf and root tissues every 4 h for one day 
(Figure 1). To confirm the microarray data and examine 
circadian rhythms of the selected clock genes in N. 
attenuata, we analyzed the transcript accumulation of 
the candidate genes in seedlings grown under different 
light conditions. N. attenuata plants were grown under 
16 h light/ 8 h dark cycle (LD) for two weeks and subse- 
quently transferred into continuous light condition (LL). 
Twenty seedlings in LD and LL were harvested every 
2 h for three days. Quantification of mRNA expression 
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Figure 1 Diurnal rhythms of putative circadian clock genes in 
Nicotiana attenuata. Expression of orthologs of the circadian clocl< 
genes were first examined using our time course microarray data 
(accession number GSE30287). Wild-type N. attenuata plants were 
harvested every 4 h for one day from leaves and roots. After RNA 
isolation, each sample was hybridized on Agilent single color 
technology arrays designed from the N. attenuata transcriptome 
(accession number GPL13527). Mean (± SB) levels of transcript 
abundance of (A) NaLHY, (B) NaTOCl, and (C) NaZTL in leaves (n = 3) 
and roots (n = 3) at each harvest time. Gray boxes depict the dark 
period of LD (16 h light/ 8 h dark). Na, Nicotiana attenuata; LHY, late 
elongated hypocotyl; TOCl, timing of cab expression 1; ZTL, 
zeitlupe. 



Yon et al. BMC Plant Biology 201 2, 12:1 72 
httpy/www.biomedcentral.com/l 471-2229/1 2/1 72 



Page 3 of 10 



was performed by quantitative real-time PGR (qPCR) 
using the gene specific primers (Additional file 1). We 
also constructed full-length of coding sequences of 
NaLHY, NaTOCl, and NaZTL (Additional file 2), based 
on the ortholog sequences available in public EST data- 
bases and Arabidopsis database. To examine evolution- 
ary relationship of circadian clock genes in plant species, 
phylogenetic trees were constructed using the UPGMA 
algorithm (Figure 2 and Additional file 3). 

Arabidopsis and poplar genomes contain two MYB 
transcription factors (LHY/CCAl and LHYl/2, respect- 
ively), which play a key role in the regulation of the en- 
dogenous clock [23,36]. However, we found only one 
oscillating AtLHY-like transcript in our cDNA library 
(Figure 2A). The transcript levels of NaLHY in LD 
peaked at Zeitgeber time 0 h (ZT 0 h) and remained low 
until the next ZT 0 h (Figures lA and 3A), consistent 
with the patterns observed in Arabidopsis and rice 
[26,37]. A/(^L//r transcript accumulation under LL main- 
tained a circadian rhythm peaking near subjective dawn 
and showed a shorter period of time with an average 
6.1 h than under LD (Figure 3A). Protein sequence of 
NaLHY shared a relatively high similarity with AtLHY 
(Identities = 42%, Positives = 54%, Gaps = 18%) and AtCCAl 
(Identities = 38%, Positives = 50%, Gaps = 22%). Phylo- 
genetic analysis indicated that NaLHY was more close- 
ly related to AtLHY than to AtCCAl (Figure 2). The 
NaLHY protein (767 aa) is larger than AtLHY (645 aa), 
AtCCAl (608 aa), and OsCCAl (719 aa) but has the con- 
served SANT domain at the N-terminal and the two 
alanine rich regions which characterize the AtLHY and 
OsCCAl proteins (Additional file 2). 

TOCl/PRRl is a member of the PRR family, which is 
composed of five oscillating genes; TOCl/PRRl, PRR3, 
PRRS, PRR7 and PRR9, Each of PRR genes has its own 
diurnal expression pattern [38]. The microarray data 
revealed that NaTOCl transcripts peaked at ZT 12 h 
(Figure IB) in accordance with the known circadian 
rhythm of the Arabidopsis TOCl [39,40]. The qPCR ana- 
lysis showed that NaTOCl transcripts peaked at ZT 
12 h under LD and the expression of NaTOCl under LL 
peaked earlier by an average of 1.3 h compared to the 
expression under LD (Figure 3B). The full-length 
NaTOCl protein sequence exhibited high similarity to 
AtTOCl/PRRl (Identities = 49%, Positives = 59%, Gaps = 
16%). Phylogenetic analysis revealed that NaTOCl is 
most closely related to AtPRRl than to other AtPRRs 
(Figure 2). The REG domain at the N-terminal and the 
CCT motif in the C-terminal of TOCl were conserved 
in N attenuata, Arabidopsis and O. sativa but the 
coiled-coil region was only found in the eudicots 
Arabidopsis and N attenuata (Additional file 2). 

The Arabidopsis genome encodes three F-box proteins 
involved in protein degradation of the clock components 
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Figure 2 Phylogenetic trees of putative circadian clocl< genes in 
N, attenuata. Phylogenetic relationships among predicted orthologs 
of (A) LHY and CCAl, (B) TOCl and PRRs, (C) ZTL and FKF1/AD03 in 
Arabidopsis thoiiono (At) and N. attenuata (Na). Full-length amino 
acid sequences were aligned and phylogenetic trees were 
reconstructed by the UPGMA method. The numbers given for each 
branch represent the numbers of amino acid substitutions per site. 
LHY, late elongated hypocotyl; CCAl, circadian clock-associated 1; 
TOCl, timing of cab expression 1; PRR, pseudo-response regulator; 
ZTL, zeitlupe; FKFl, flavin-binding kelch repeat F-box 1; ADO, adagio. 



[18,19,21]. We found two ZTL orthologous genes 
(NaZTL and NaFKFllADOS) in our cDNA library and 
named them according to the phylogenetic analysis 
(Figure 2). NaZTL transcripts under LL and LD were 
largely unchanged (Figures IC and 3C), consistent with 
the results from Arabidopsis [41]. In contrast, NaFKFI 
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Figure 3 Circadian rhythm of the clock gene expression in 
N. attenuata. Mean (± SE) levels of relative transcript abundance of 
(A) NaLHY, (B) NaTOCl, and (C) NaZTL in wild-type N. attenuata 
seedling (n = 3) at each harvest time for two days under long day 
condition (LD, gray lines, 16 h light/ 8 h dark) and continuous light 
condition (LL, black line). Gray boxes depict the dark period of LD. 
A N. attenuata elongation factor gene was used as a control for 
constitute expression. 



transcripts showed a clear circadian rhythm with peaks 
at ZT 12 h under LD (Additional file 4). The period of 
NaFKFl expression in LL was shortened by an average of 
2 h (Additional file 4). NaZTL had protein sequence 
similarity to AtZTL (Identities = 79%, Positives = 86%, 
Gaps = 5%) and NaFKFl with that of AtFKFl (Identities = 
78%, Positives = 85%, Gaps = 4%). NaZTL and NaFKFl 
proteins contained LOV/PAS and F-box domains at 
their N-terminals and Kelch repeats region in their 
C-terminals as shown in orthologs in Arabidopsis 
(Additional file 2). 

Ectopic expression of N. attenuata clock genes in 
arabidopsis 

To examine the functional conservation of the clock com- 
ponents in N, attenuata, we produced Arabidopsis (Col-0) 



transgenic plants ectopically expressing NaLHY and 
NaZTL genes under the control of Cauliflower Mosaic 
Virus (CaMV) 35S promoter {35SNaLHY and 35S: 
NaZTL), The circadian clock entrains hypocotyl ceU 
elongation to light-dark cycles, and thus plants overex- 
pressing either Arabidopsis LHY or ZTL have elongated 
hypocotyls compared to wild-type plants [37,42]. We mea- 
sured the hypocotyl lengths of 35SNaLHY (Figure 4A) 
and 35SNaZTL transgenic plants (Figure 4B). Two inde- 
pendent 35SNaLHY lines showed a pronounced increase 
in hypocotyl length compared to wild-type, Col-0, as did 
the two independent 35SNaZTL lines (P<0.05, one-way 
ANOVA with Fishers post hoc test). 

The ability to perceive seasonal changes by the circa- 
dian clock is required for the successful transition from 
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Figure 4 Ectopic expression of NaLHY and NaZTL in A. thaliana. 

Arabidopsis transgenic plants harboring 35S\NaLHY and 35S\NaZTL 
constructs are sliown at seedling (A, B) and flowering (C) stages. 
Mean (± SE) values of hypocotyl lengths in wild-type Col-0 (n = 1 1), 
35S:NaLHY-2 (n = 1 2), 35S:NaLHY-3 (n = 1 1 ), 35S:NaZTL-] (n = 1 3), and 
35S:NaZTL-2 (n = 13). Different letters (a and b) reflect significant 
differences among the lines (P<0.05, one-way ANOVA with Fisher's 
post hoc test). (C) The percentage of flowering plants and number of 
rosette leaves of Col-0 (n = 1 6), 35S:NaLHY-2 (n = 32), 35S:NaLHY-3 
(n = 31), and 35S:NaZTL-] (n = 32), 35S:NaZTL-2 (n = 30) when 
inflorescence elongation started. 
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vegetative to reproductive stages. Knocking-out a mor- 
ning element, Arabidopsis LHY or CCAl results in an 
early-flowering phenotype and, in contrast, their overex- 
pressing lines result in a late-flowering phenotype in 
Arabidopsis [12,36,37]. The overexpression oi Arabidopsis 
ZTL also results in a late-flowering phenotype under LD 
[42]. To examine the eflect of NaLHY and NaZTL on 
flowering time regulation in Arabidopsis^ we scored the 
rosette leaf number at a time when initial flowering was 
observed in plants grown under LD. Both 35S:NaLHY 
and 35S:NaZTL plants exhibited delayed flowering phe- 
notypes and increased rosette leaf numbers (Figure 4C). 

Interaction between NaTOCI and NaZTL protein 

F-box protein ZTL plays a key role in protein turnover of 
the clock components, including TOCl by direct protein- 
protein interactions [18]. We determined whether NaTOCI 
interacts with NaZTL by yeast two-hybrid analysis. In 
addition, we tested inter-species interactions of AtTOCl- 
NaZTL and AtZTL-NaTOCl. Full-length NaTOCI pro- 
teins (AtTOCl) fused with the GAL4 activation domain 
(AD), and NaZTL (AtZTL) fused with the GAL4 binding 
domain (BD) were co-expressed in yeast cells (see 
Methods). The yeast cells expressing both NaTOCI and 



NaZTL constructs grew well in the absence of leucine, 
tryptophan, histidine, and adenine, showing the positive 
interaction, consistent with Arabidopsis TOCl and ZTL 
interaction (Figure 5A). The transformants expressing both 
NaTOCI and NaZTL also exhibited strong |3-galactosidase 
activity compared to control yeast cells carrying empty vec- 
tors or only one of the constructs (Figure 5B). Inter- 
species protein interactions of AtTOCl -NaZTL and 
AtZTL-NaTOCl were as strong as the intra-species inter- 
actions of TOCl and ZTL (Figure 5). 

Silencing of NaTOCI in N. attenuata 

To examine the conserved and unique functions of cir- 
cadian clock genes in N. attenuata, we first silenced 
NaTOCI expression by constitutive overexpression of 
A/<^ rOCi -specific inverted repeat (ir) sequences (irTOCi 
lines) [43,44]. TOCl in Arabidopsis plays a key role in 
flowering time regulation. The semi-dominant tocl-1 
mutant displays a late-flowering phenotype under LD 
and an early-flowering phenotype under short day con- 
ditions (SD) [39] and the tocl-2 loss-of- function mutant 
shows no significant difference of flowering time under 
LD but an early-flowering phenotype under SD [45]. Dur- 
ing the screening of the silenced lines, we clearly observed 
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Figure 5 Protein interactions between TOCl and ZTL from N. attenuata and A. thaliana using a yeast two-hybrid assay. (A) Growth of 
the yeast cells carrying prey and bait constructs indicated on the top of each panel. The ZTL proteins were fused to the GAL4 DNA binding 
domain (bait), and the TOCl proteins were fused to the GAL4 activation domain (prey). The yeast cells can grow on QD medium when bait and 
prey proteins physically bind. (B) Mean (± SE) levels of |3-galactosidase activity of the yeast carrying prey and bait constructs. BD, empty vector 
expressing binding domain; AD, empty vector expressing activation domain; -LW, synthetic dropout (SD) yeast growth medium lacking leucine 
and tryptophan; QD, SD medium lacking Leu, Trp, histidine, and adenine. 
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the late-flowering phenotypes of several independent 
ivTOCl lines under LD (Figure 6). Two independent 
lines were used to measure flowering time and rosette 
leaf number when the first flower opened (Figure 6C). 
The silencing efficiency of independent lines harboring 
the RNAi construct correlated strongly with the flowe- 
ring phenotype (Figure 6B). In one homozygous ivTOCl- 
162 line, screened on the basis of antibiotic resistance, 
but in which the TOCl gene expression was not silenced, 
presumably due to methylation of the ir part of the 
T-DNA or an unknown insertion effect, displayed a WT 
flowering pattern (Figure 6). 

Discussion 

Since the clock genes of Arabidopsis have been identi- 
fied, the clock mechanism of Arabidopsis has been 
extended to other dicotyledonous plants; soybean [24], 
chestnut [29], Brassica rapa [46], poplar [23] and also to 
monocotyledonous plants; rice [25], maize [22], duck- 
weed [28]. Here, we identified three core circadian clock 



genes in a wild tobacco, N. attenuata. Analyses of the 
circadian rhythms of transcript accumulation and protein 
similarity allowed us to identify orthologs of LHY, TOCl, 
and ZTL of Arabidopsis in N. attenuata. Protein interac- 
tions of TOCl and ZTL in Arabidopsis were also con- 
served in N. attenuata. In addition, ectopic expression of 
NaLHY and NaZTL in Arabidopsis confirmed the func- 
tional conservation of LHY and ZTL in N, attenuata. We 
also demonstrated that NaTOCl in N, attenuata plays an 
important role in the regulation of flowering time. 

The number of LHY (or CCAl) orthologous genes and 
the timing of gene duplication events differ among plant 
species [23]. One common ancestor of LHY/ CCAl inde- 
pendently duplicated in monocots and eudicots. Populus 
nigra and P, trichocarpa contain two LHY orthologs and 
the monocots, rice and Sorghum bicolor contain one 
CCAi -like gene in their genome [23,47]. Gene duplication 
events of LHY I CCAl in popular and Arabidopsis would 
be expected in N, attenuata but we were able to find only 
one LHYICCAl like gene in N. attenuata. However, the 
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Figure 6 Late-flowering piienotypes of NaTOCl silenced N. attenuata. (A) Phenotypes of transgenic plants silenced in NoTOCl gene 
expression (ir70C/) with an inverted repeat (ir) RNAi construct. (B) Mean (± SE) levels of NoTOCl (n = 3) transcript accumulation at zeitgeber time 
(ZT) 12 in wild-type and in several ir70C/ lines. (C) Time on weeks ofWT and 2 selected irTOCl lines to reach flowering stage (n = 5). (D) Mean 
(± SE) number of leaves at flowering stage (n = 5). 
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second ortholog could be missing in our current cDNA 
library and we plan to perform deep sequencing of the 
transcriptome and microarray analysis of plants grown 
under a variety of conditions in the future to clarify the 
evolutionary relationship of N. attenuatds LHY orthologs. 

While components of the endogenous clock and their 
associated circadian clock mechanisms have clearly been 
maintained across diverse plant species, clock-mediated 
signaling has evolved in response to differing selection 
pressures. A perennial plant, chestnut {Castanea sativa) 
has orthologs of AtLHY and AtTOCl in its genome and 
the pattern of transcript accumulation of CsLHY and 
CsTOCl is similar to that of Arabidopsis in LD at 22°C 
[29]. However, in winter condition, both CsLHY and 
CsTOCl transcripts lose their diurnal rhythms and 
maintain high levels of transcripts, which may be asso- 
ciated with the induction of winter dormancy [29] . Even 
within a single species. A, thaliana, genetic variation in 
the clock components plays a critical role in adapting 
ecotypes to their local environment [48]. Genetic 
variation in the PRR genes is associated with local adap- 
tation seen in the differential expression of quantitative 
trait loci in 150 Arabidopsis ecotypes [48]. In addition, 
the degree to which TOCl regulates flowering time dif- 
fers among Arabidopsis ecotypes. A semi-dominant 
tocl-1 mutant of the Arabidopsis C24 ecotype displays 
late a flowering phenotype in LD, whereas the same 
TOCl mutation in the Landsberg ecotype results in no 
change in flowering time compared with WT plants 
under LD [39]. The toe 1-2 mutant shows also no change 
in flowering time under LD [45]. However, silencing 
TOCl in N, attenuata confers a late-flowering pheno- 
type under LD, which may be due to the longer life span 
of N, attenuata (about 3 months) or different circadian 
clock functions have evolved under various environmen- 
tal pressures. In future research, we plan to measure the 
flowering time of ir TOCl lines under SD to examine the 
light sensitivity of this transgenic line. 

We have investigated the ecology of N. attenuata in 
its native habitat for more than twenty years. During this 
period, we have observed interesting diurnal rhythmic 
traits of N. attenuata and time-of-day dependent eco- 
logical interactions. For example, N. attenuata interacts 
with different groups of herbivores which are either day- 
active (such as grasshoppers, mirids and Manduca 
larvae) or night-active (such as noctuid larvae and tree 
crickets) and produces different chemicals that function 
as direct defenses against these herbivores or function as 
indirect defenses and attract of predators of the herbi- 
vores [49]. Recently, we showed that tissue specific diur- 
nal rhythm of metabolites and its related transcripts in 
N, attenuata changes in response to herbivore attack of 
a specialist, M. sexta larvae [35]. More than 15% of total 
metabolites that we measured in leaf and root shows 



diurnal patterns and some of them have been demon- 
strated to function as plant defenses against herbivore 
attack. Goodspeed et al, [50] have recently reported that 
feeding behavior of Tricoplusia ni is predicted by the 
circadian clock in its host plant Arabidopsis and it in- 
creases anti-herbivore defense of Arabidopsis, All of 
these interactions provide a rich arena in which to ex- 
plore the molecular mechanism of how the circadian 
clock regulates plant-insect interactions. 

Conclusions 

We identified three core circadian clock components in 
N, attenuata based on the gene expression and pheno- 
typic alterations in lines silenced or overexpressed in the 
components. This work provides the foundation for the 
manipulation the ecological roles of the circadian clock 
in nature. As jet travel has revealed the depth of circadian- 
regulated processes in humans, circadian mutants of 
N, attenuata will be used to unravel the ecological func- 
tions of the clock in plant-environment, plant-plant, and 
plant-insect interactions in nature. 

Methods 

Plant growth condition 

For all experiments we used Nicotiana attenuata Torr. 
Ex. Wats (Solanaceae) plants (31^^ inbred generation), 
wild-type (WT) originating from a population in Utah. 
Seeds were sterilized and germinated on Petri Dishes 
with Gamborgs B5 media as described in [43]. Petri 
dishes with 20 seeds were kept under long day condi- 
tions (LD, 16 h light/ 8 h dark). 

To examine the free running period, one group of 
seedlings 15 days after germination was moved into con- 
stant light conditions for the 3 days of sampling and the 
other group of seedlings was grown under LD. Three 
biological replicates (10 seedlings pooled per replicate) 
were harvested every 2 h for 3 days and immediately fro- 
zen in liquid nitrogen. 

RNA isolation and gene expression 

Total RNA was extracted using TRIZOL reagent and 
cDNA was synthesized from 500 ng of total RNA using 
RT- minus kit (Fermentas, Burlington, Canada). The quan- 
titative RT-PCR analyses were performed on a Stratagene 
MX3005P (Agilent Technologies, Santa Clara, CA, USA) 
employing SYBR Green kits (Eurogentec, Cologne, 
Germany). Primers were designed based on sequences 
from Arabidopsis thaliana retrieved from TAIR website 
and a N, attenuata transcript library (NCBI GEO Data- 
base accession number GSE30287). The sequences of 
qRT-PCR primer pairs for NaLHY, NaTOCl, NaZTL and 
NaFKFl are listed in Additional file 1. Transcript abun- 
dance expressed relative to the expression of N, attenuata 
ELONGATION FACTOR gene. 
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Phylogenetic analysis 

The identity of the circadian genes NaLHY (NCBI acces- 
sion number JQ424913), NaTOCl (Accession number 
JQ424914), NaZTL (Accession number JQ424912) and 
NaFKFllADOS (Accession number JQ424915) was deter- 
mined by sequencing. A standard PGR was performed to 
obtain the full length cDNA and the list of primers used 
for this analysis is in Additional file 1. PGR products were 
subcloned for amplification using a GloneJET PGR 
Gloning kit (Fermentas). 

The amino acid sequences of the circadian genes 
NaLHY, NaTOCl, NaZTL and NaFKFllADOS were 
deduced from cDNA sequences and aligned using the 
Geneious program V5.3 (http://www.geneiousxom). The 
numbers of amino acid substitutions were estimated 
by a Jukes-Gantor model using a BLOSUM 62 matrix, 
through a global alignment with free end gaps option. A 
phylogenetic tree was reconstructed by the Unweighted 
Pair Group Method with Arithmetic Mean (UPGMA) 
method. These analyses were performed using Geneious 
software V5.3. 

Overexpression of NaLHY and NaZTL In arabidopsis 

All A. thaliana lines used were in the Gol-0 background. 
Plants were grown in a controlled culture room at 22°G 
with a relative humidity of 55% under long-day (LD) 
conditions (16-h light/8-h dark) with white light illumin- 
ation (120 (imol photons/m^s) provided by fluorescent 
FLR40D/A tubes (Osram). To produce transgenic plants 
overexpressing the NaLHY 2ind NaZTL genes, full-length 
cDNAs were subcloned into the binary pB2GW7 vector 
under the control of the GaMV 35S promoter (Invitro- 
gen, Garlsbad, GA, USA). Agrobacterium tumefaciens- 
mediated Arabidopsis transformation was performed 
according to a modified floral dip method (Glough and 
Bent, 1998). T2 transgenic plants harboring a single 
T-DNA insertion were used in subsequent assays includ- 
ing hypocotyl length and flowering time measurements. 
Transformation with NaTOCl was carried out but for 
unknown technical reasons we were not possible to re- 
generate plants from it. Selected lines were checked by 
RT-PGR for the NaLHY and NaZTL overexpression 
(Additional file 5), using A, thaliana tubulin (TUB) as a 
reference gene. 

Yeast-two-hybrld 

Yeast two-hybrid assays were performed using the 
Matchmaker™ system (Glontech, Palo Alto, GA, USA). 
Primers used for amplif)^ing NaTOCl, AtTOCl, NaZTL, 
and AtZTL were described in Additional file 1. Each 
RT-PGR product was digested with restriction enzymes 
{Ecom, Xmal for NaTOGl and AtTOGl, Ncol, Xmal for 
NaZTL, Ndel, BamHl for AtZTL). Digested full-length 
transcripts of TOCls and ZTLs were subcloned into 



pGADT7 and pGBKT7, respectively. The yeast strain 
AH109 (leucine-, tryptophan-, histidine-, adenine-) con- 
tained lacZ reporter gene was co-transformed with the 
indicated vectors in Figure 5A. Transformation was con- 
ducted according to the manufacturers instructions 
(Glontech). Single colonies obtained on growth medium 
lacking Leu, Trp were inoculated on a medium without 
Leu, Trp, His, Ade and used to measure p-galactosidase 
activity described in the instructions (Glontech). 

Silencing of NaTOCl In N. attenuata 

A sequence fragment of NaTOCl cDNA was inserted 
into the pSOLS transformation vector as an inverted re- 
peat construct driven by the GaMV 35 S promoter 
[43,44]. The NaTOCl vector was transformed into 
N attenuata WT plants using Agrobacterium tumefa- 
ciens mediated transformation; ploidy was determined 
on To plants as described by Gase et al. [44], allowing 
for the selection of only diploid transformed lines. 
Homozygosity was confirmed on T2 plants by hygromy- 
cin resistance screening after which 6 transformed lines 
were selected and transferred to the glasshouse for fur- 
ther growth at conditions described in Kriigel et al. [43]. 
Gene expression levels of NaTOCl were determined by 
qPGR from tissue of selected T2 plants and wild- type 
plants collected at ZT 12 h. 

Additional files 



Additional file 1: List of primers used for transcript profiling and 
full-length cloning of circadian clock genes in N. attenuata and A. 
thaliana. 

Additional file 2: Protein alignments of circadian clock genes in N. 

attenuata, Arabidopsis and rice. Full-length amino acid sequences were 
aligned using the Geneious software. 

Additional file 3: Phylogenetic trees of circadian clock genes in 
several plant species. Phylogenetic trees of (A) LHY/CCAl, (B) TOCl, (C) 
ZTL and (D) FKF1/AD03. Full-length amino acid sequences were aligned 
using the Geneious software. Phylogenetic trees were constructed with 
available sequences from 3 major taxonomical groups: eudicots, 
monocots and one moss Selaginello moellendorffii for TOCl and ZTL 
trees. Unweighted Pair Group Method with the Arithmetic Mean 
(UPGMA) method from the numbers of amino acid substitutions by 
applying the Jukes-Cantor model was used. The scale bar represents the 
number of amino acid substitutions per site. Ac, Allium cepa; At, 
Arabidopsis thaliana; Cs, Castanea saliva; Gm, Glycine max; In, Ipomoea nil; 
Mc, Mesembryanthemum crystallinum; Na, Nicotiana attenuata; Os, Oryza 
sativa; Pn, Populus nigra; Pt, Populus trichocarpa; Pv, Phaseolus vulgaris; Sb, 
Sorghum bicolor; SI, Solanum lycopersicum; Sm, Selaginella moellendorffii; 
Vv, Vitis vinifera; Ta, Triticum aestivum; Zm, Zea mays. 

Additional file 4: Circadian rhythm of the NaFKFl/ADOS expression 
in N. attenuata. (A) Mean (± SE) levels of transcript abundance of 
NaFKFl/AD03 were examined using our time course microarray data 
(accession number GSE30287). Wild-type N. attenuata plants were 
harvested every 4 h for one day from leaves (n = 3) and roots (n = 3). 
After RNA isolation, each sample was hybridized on Agilent single color 
technology arrays designed from the N. attenuata transcriptome 
(accession number GPL13527). (B) Mean (± SE) levels of relative transcript 
abundance of NaFKFl/AD03 in wild-type N. attenuata seedling (n = 20) at 
each harvest time for two days under long day condition (LD, gray lines. 
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16 h light/ 8 h dark) and continuous light condition (LL, black line). 
Transcript levels were determined by qPCR and N. attenuata Elongation 
Factor was used as reference gene. Gray boxes depict the dark period 
of LD. 

Additional file 5: Overexpression of NaLHY and NaZTL transcripts in 

Arabidopsis transformants. Transcript levels were measured by RT-PCR. 
Primers (see Additional file 1) were designed for the specific detection of 
NaLHY or NaZTL transcripts, not for Arabidopsis LtiY or ZTL. Plants grown 
under LD were harvested at ZT 8 h. A tubulin gene (TUB) was included 
as a control. 
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